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I ■ Abstract 

^: 

r-| [ We discuss two ways of deriving the fluctuation-dissipation theorem (FDT) for soft fermion 

', excitations in a hot non-Abehan plasma being in a thermal equilibrium. The first of them is 

^ ■ based on the extended (pseudo) classical model in describing a quark-gluon plasma suggested in 

; [YuA. Markov, MA. Markova, Nucl. Phys. A 784 (2007) 443], while the second one rests on 

2 , the standard technique of calculation of the FDT for thermodynamically equilibrium systems. 

G^ ' We show that full accounting all subtleties that are common to the fermion system under 
O 

Q\ consideration, results in perfect coincidence of thus obtained FDTs. This provides a rather 

o 



X 



strong argument for the validity of the pseudoclassical model suggested. 
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1 Introduction 



In our two previous papers |TP] we have proposed an approach to construction of the generahzed 
(pseudo) classical theory for a unified description of a broad spectrum of interaction processes 
of soft and hard quark-gluon plasma excitations obeying both Fermi and Bose statistics. Set- 
ting into consideration Grassmann (anticommuting) color charges 9 = (6'*) and O'' = (^^^*), i = 
1, . . . , Nc, along with usual (commuting) classical color charge Q = (Q"), a = 1, . . . , — 1, of 
hard particle, is a principle feature of this extended theory. A number of examples of an ex- 
plicit constructing the probabilities for various scattering processes involving the bosonic and 
fermionic hard and soft modes, has been given. Also possible applications of the theory we 
have developed to an important problem of the calculation of energy losses of fast Yang-Mills 
particles passing through a hot QCD matter, were considered. 

However, in spite of the reasonableness of presuppositions in the construction of the gen- 
eralized theory mentioned above and also obtaining rather nontrivial results within the limits 
of this theory, we would like to have independent, physically transparent and simultaneously 
a simple way for critical examination of the approach suggested in [HE]- This independent 
testing must be founded on some fundamental, i.e. model- independent, physical principle. As 
such a principle in the present work we have chosen the fluctuation- dissipation theorem which 
as is generally known [3j, admits an abstract formulation convenient for applications of this 
theorem to fields of equilibrium thermal fluctuations of any physical nature: electromagnetic, 
mechanical, temperature, entropy and so on. The FDT already for a long time, is an integral 
part in courses on the statistical mechanics |1H9] (see, also recent review |10)). 

A distinctive characteristic of the problem under consideration, is the fact that an external 
perturbation described by a Hamiltonian H^, changes the number of particles and antiparticles 
in the system, i.e., in other words, an operator of the total number of particles N (more precisely, 
the difference between the number of particles and antiparticles) does not commute with the 

(but it commutes with a Hamiltonian of the many-particle system Hq). And though such a 
case has already been discussed in the general statement in textbook materials (see, e.g., [6l|9]), 
nevertheless as far as we know, any specific physical situation, where this circumstance would 
play crucial role, has not been considered. As a consequence of this fact, a concrete expression 
for the fluctuation-dissipation theorem in which this noncommutativity could be manifested by 
obvious fashion, has not been given anywhere. For this reason, in the present paper we have 
found pertinent to give the comprehensive derivation of the FDT relations for the physical 
problem at hand. 

The paper is organized as follows. In Section 2, by using a simple classical model of the 
quark-gluon plasma, we derive the spectral density for thermal bosonic fluctuations. We show 
that the expression obtained exactly coincides with that following from the corresponding 
fluctuation-dissipation theorem. Section 3 is concerned with deriving relevant spectral den- 
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sity for thermal fluctuations of the quark- ant iquark field in the medium within the extended 
(pseudo) classical model of the QGP. Sections 4 and 5 are devoted to determining the FDT rela- 
tions for the volume field density of the generalized forces and the stationary and homogeneous 
soft fields of the system in question, obeying Fermi-Dirac statistics. In the former section we 
obtain the FDT for fluctuations of the hard sources r^^ and f/^, while in the latter a similar the- 
orem is derived for fluctuations of the soft spinor fields ipl^ and ipl^. At the end of Section 5 we 
discuss the condition of matching two expressions for the soft-quark spectral density obtained 
within the pseudoclassical model of QGP and the fluctuation-dissipation theorem. 



2 Correlation function for soft Bose-excitations 

Before we proceed to calculation of the spectral density for fermionic fluctuations and to deriva- 
tion of the relevant fluctuation-dissipation theorem, we would like to consider briefly similar 
questions for much-studied case of fluctuations of a gauge field in the medium. There are 
no new results to be reported here. The main goal of this section is to outline the principle 
points of obtaining the FDT for soft bosonic fluctuations. Further, in subsequent sections 
we shall follow these points in deriving a similar FDT for the fermion degree of freedom of a 
non-Abelian plasma. A somewhat different approach to obtaining the fluctuation-dissipation 
theorem for bosonic sector of plasma fluctuations can be found in the papers by D.F. Litim 
and C. Manuel HI]. 

As a first step we shall calculate the spectral density for bosonic fluctuations in the hot 
non-Abelian plasma within the classical model suggested by U. Heinz more than 20 years 
ago [T2HI1]. The second step is a comparison of the obtained expression for the spectral den- 
sity with a similar expression resulting from the FDT for fluctuations of a gauge field. 

According to jT3], the soft gluon field induced by a hard test particle (which is centered 
at the position xg) in the momentum representation is 

A^(fc) = - *V^,{k)r{k- Xo), (a;, k) (2.1) 

where in turn the current of the hard color-charged particle reads 

r{k-^) = j^^v''Q''5{vk)e-^'^--\ ^ = (l,v). (2.2) 

Here, *V^^{k) is the gluon propagator in the hard thermal loop (HTL) approximation^; v and 
are the velocity and classical color charge of the hard particle, respectively. In the simplest 

^ In the high-temperature QCD within the hard thermal loop approximation it is customary to denote all the 
effective (i.e. HTL-resummed) propagators and vertices by *!?, *iS', *r and so on, with the aim to distinguish 
them from the exact quantities I?, 5, F, . . . (see, for example, the review [15J). We shall maintain these notions 
throughout our paper. As usually we put the sign of complex conjugation " * " to the right of a quantity. 
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approximation the last two quantities are taken time-independent ones. 

Further, we consider the correlation function (^A*^ (k) A'l{k')'^ with gauge field (12. ip . We 
need first to define what is meant by the angular brackets of averaging. The average procedure 
involves the integration over all possible initial positions xq of the hard particle. Then, the 
average involves the integration with respect to color charge with the measure dQ. An explicit 
form of the latter is written out in Appendix A. Besides, the average over the distributions of 
hard particles: quarks (Q), antiquarks (Q) and gluons (G) being in thermal equilibrium, should 
be added. Thus we can write 



X {*V^^,{k)r {*V,.'{k')) l^v^'v'''6{v ■ k)6{v ■ k'). 



(2.3) 

■v^'^/Mv ■ k)^(v ■ k'). 

47r 

The overall factor 2 takes into account that hard gluons and massless hard (anti) quarks have 
two helicity states. The integration over Xq is trivially performed 

j rfxo e^C^-i^')- ^0 = {2nf6{k - k') . (2.4) 

Furthermore, considering the last expression, in the second line of (12. 3p . we have 

6{v ■ k)6{v ■ k') = 6{u - uj')6{u - V • k). (2.5) 

The integral over the color charge is defined by a formula 

J dQQ^Q^= (^j^ dA = Nl - 1. (2.6) 

The exact definitions of the constants C^"^ in terms of the group invariants, are given in Ap- 
pendix A. It takes a little more efforts to analyze the integral over a solid angle. Let us 
present it in the form of an expansion in terms of the basis of Lorentz-covariant matrices fixed 
in the rest frame of the medium 

' dVL^ 



^ v^v^'S^v ■ k) = atP^'{k) + ai Q^'{k) + . . . , 
47r 

where at^i are unknown for a while constants; P^*^ and Q^^ are transverse and longitudinal 
projectors whose explicit forms are given in Appendix B. The dots designates contributions of 
the other possible tensor structures. 

Let us introduce the spectral density for soft bosonic fiuctuations through a relatione 

{Al\k)Al{k')) = {A*;AlU5{k - k'). (2.7) 



On the right-hand side of the definition of the spectral density (|2.7p . the factor (27r)^ is absent since we use 
the following formulae for the direct and inverse Fourier transforms: 

F{x) = f F{k)e-'''-''dk, F{k) = f F{x)e'''-'- 



i.e., the factor l/(27r) attaches to the inverse transform. 
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Substituting Eqs. (12.41) - (12.61) into (12.31) . taking into account a decomposition of the gluon prop- 
agator in the covariant gauge, Eq. (B.l), and the properties of projection operators (B.4), we 
find a required form of the spectral density within the hmits of the classical model 

^(0- 



2(4^"^E (-) 



Y^a /\ b 

C=Q,Q,G 



x\atP^,{k)\*A\k)\'+ aiQ^,{k)\*A\k)\ 



On the other hand, the spectrum of soft gluon modes in an equilibrium QGP can be defined 
by means of the fluctuation-dissipation theorem. Here, we can use the well-known expression for 
the spectral density of fluctuations of an electromagnetic field [T6|IT7] with a minimal extension 
to the color degrees of freedom 

{A*;KU = -jiy.^^'' (^^^){^-^^) - (2-9) 

where 

is the mean energy of a quantum bosonic oscillator and kB is the Boltzmann constant. As 
for the factor 1/(27?)'^ on the right-hand side of (12. 9p . see footnote [2j In the semiclassical 
approximation the gluon (retarded) propagator should be taken in the HTL-approximation, 
Eq. (B.l). Then in the limit h ^ 0, the expression (12.91) results in 

= ji^J^'(^){P,^(k)lmCA-'\k)r\*A\ 

In the rightmost-hand side of the previous equation we have taken into account the obvious 
identity 

Im*A*''(A;) = Im(*A-i*''(A;))* |*A*''(A;)| 2. 
The requirement of consistency of (12. 8p with expression (12. lip leads to the conditions 



9^ 
2 ^ 



(27r) 



C=Q,Q,G 



dA .r 27r2 



For definiteness we will consider the condition of consistency for the longitudinal part of the 
expressions (12. 8p and (12. lip . By virtue of (B.l) and (B.2), we have 



,2 



Im(*A-^'(u;, k))* = uTfOik" - 

2 kr 



where m|, is the Debye screening mass. 
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Furthermore, in the covariant gauge the longitudinal projector Q^'^{q) in the rest frame of 
the heat bath, is equal to 

and therefore, a; is defined as follows: 

Here, we have used the integration formulae over a solid angle (B.IO). Based on the above- 
mentioned, the condition of consistency fl2.12p will be identically fulfilled for the longitudinal 
part (and as it is easy to see, it is also valid for the transverse part), if the Debye screening 
mass is defined by a relation 



2 C=Q,Q,G\ "A /J ^TT 

= ■ (2-13) 

It should be noted, however, that the expression on the right-hand side of f l2.13p represents 
a definition for the classical Debye mass squared, while on the left-hand side m^, stands for 
the quantum Debye mass squared. Let us recall that the latter appears in the definition of 
the HTL-resummed gluon propagator (B.1)-(B.3). By virtue of the fact that the classical 
Debye mass differs from the quantum one, this circumstance results in violation of the FDT 
under this 'naive' approach. This fact was first noted in the works [11]. The matter is that 
in the average ( 12. 3 p we use the procedure of the classical statistical average, usually accepted 
in a high-temperature Abelian plasma. However, for a hot non-Abelian plasma it is not quite 
true. To derive the correct quantum value for the Debye mass (and thereby to recover the 
FDT) following [llj, we should replace the classical average in (12. 3p by the quantum one that 
holds for the systems close to thermal equilibrium. In this case it reduces to the following 
replacements in (12.31) 



where Nf is the number of active (massless) quark fiavors. Then taking into account the last 
circumstance we will have the following expression for the Debye mass instead of ( I2.13P 



.(Q) 



+ 



(2.15) 



For values of the quadratic Casimirs C2 specified in Appendix A, the expression fl2.15p correctly 
reproduces the quantum Debye mass squared as it was defined in the high-temperature QCD 

m- 

Thus the simple model for the classical description of the hot non-Abelian plasma with the 
known proviso enables us to exactly reproduce the FDT relation for soft bosonic fluctuations. 



3 Correlation function of soft Fermi excitations 

We proceed now to the problem of calculation of the correlation function for thermal fermionic 
excitations. Within the model suggested in [lj,i2j for a description of the soft and hard fermion 
degrees of freedom of the system under consideration, the soft spinor fields ip and induced 
by a free spin-1/2 hard test particle, are 

'^ai-q) = V'pi-i; Xo) *Sl3a{-q), 

where a color source r/^ of the particle has the form 

Vliq; xo) = ^ 6{v ■ q) e'^^-\ (3.2) 

and f]^{—q; xq) = ri^^{q; Xo)7''. Here, 6^ is a Grassmann color charge and X/3 is a c-number spinor 
describing a polarization state of the spinning particle. 

Let us consider the correlation function Q''))- First of all we will discuss a 

question what is meant by the angular brackets of averaging. As in the bosonic case the average 
includes the integration over every possible initial positions xq of the hard test parton. Besides, 
the average also includes the integration over Grassmann color charges with the measure dOdO^ 
specified in Appendix C, Eqs. (C.2) and (C.5). Finally, it is necessary to define the correct 
statistical factor. As was shown in yj, the relevant statistical factor here must be 

E./p^[4"+/r]^. (3-3) 

C — Qi Q 

From the above-mentioned and the formula of integration with respect to Grassmann color 
charges (C.6), we obtain 



(27r)6 ^ I 



27r2 ^ ^ 

(3.4) 



We notice that unlike the bosonic case now it is impossible to write down the integral over the 
magnitude of hard momentum |p| in the form of an independent multiplier since the spinors x 
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and X in the second line of the above expression also implicitly depend on |p|, see below. 

Let us analyze the integrand in (13.41) with the quark propagators. We shall restrict ourselves 
to the case of fully unpolarized state of (massless) hard test particle. As was shown in Appendix 
C of the paper jjj, the following replacement 



is valid. Then we get 



*S{q)x)ArS{-q)), = ^^CS{q){v^yS{-q))^,. (3.5) 



An explicit form of the quark propagator in the HTL-approximation is given in Appendix B, 
Eq. (B.5). Within this approximation the propagator represents an expansion in terms of the 
spinor projectors /i±(q), whose explicit forms are defined by equation (B.6). Our first task is 
to present the matrix (t> ■ 7) = 7° — v • 7 entering into the right-hand side of (13. 5p in the form 
of a similar expansion. For the 7*^ matrix we have 

7° = /^+(q) + /^-(q)- 

Then we write a scalar product v • 7 = v'^S^^'j^ in the form of an expansion in three-dimensional 
transverse and longitudinal projectors with respect to the vector of momentum q 

^^(5^^'-gY)7^ + (v-q)(7-q), q = rj- 

|q| 

Here, in the second term we can also write 

(7-q) = -(/^+(q)-/^-(q))- 

From what has been said, we finally derive 

(t; ■ 7) = (1 + V • q)/i+(q) + (1 - v ■ q)/i_(q) - 7 ■ (q x (v x q)). 

The last term on the right-hand side of the above expression upon integrating over a solid 
angle dQ^^/ATC in view of the function 6{u — v ■ q), will turn into zero. Therefore, subsequently 
this term will be omitted. Substituting the expression for [v ■ 7) obtained by this means into 
(13. 5p . taking into account the structure of the quark propagator and the properties of the h±{q) 
matrices (B.9), we derive instead of (13. 5p 



{*S{q)x)®{x*S{-q)) 

1 

The angular integration dQ^ in (13. 4p is easily performed employing Eq. (B.IO). We present 



(1 - V q)/.+ (q) rA+(g)|2 + (l + v q)/._(q) \*A4q)\' 



the final expression for the correlator (13. 4p in the form similar to (12. 7p . i.e., 
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where 

mu-^''^.^\ (3-6) 



l-^)(/^+(q)WrA+(g)r + (^l + ^)(/._(q)WrA_(g)P 



is the spectral density for soft quark fluctuations, and 

rf|p|. (3.7) 



Wo 




AC) , AG) 
/|p| +/|p| 



From the last expression we see that the quantity Uq in exact coincides with the fermion plasma 
frequency squared if we set 

^(Q) _ ^(Q) ^ _ ^3_g^ 

Recall that the constant Cf""^ is connected with an integral of motion for a system of dynamical 
equations (C.l) and it is the only free parameter in the model under consideration. The above- 
mentioned requirement uniquely fixes this constant. It is interesting to note that with the 
choice of (13. 8 P the constant C^P correct to a sign, coincides with the quark quadratic Casimirs 
(A. 3) by virtue of the identity 

CFN, = TFdA. 

Thus this choice looks quite reasonable. 

In addition, one further remark is in order. In contrast to the expression for the Debye 
screening mass (I2.13p . we have immediately derived the expression for the fermion plasma fre- 
quencjlf] (13. 7p . as it appears in calculating within the high-temperature QCD [T9|[20]. This 
circumstance suggests that (13. 3p represents not classical, but quantum average. Strictly speak- 
ing, the statistical average should be written in the following form: 



C=Q,Q 



/p'^[/f(l-/»)+/?'(l + /fO](lF. (3.9) 



which is 'fermionic' analog of the quantum average (I2.14p . A somewhat unusual mixed form of 
arrangement of the distribution functions for thermal particles is caused by the following fact: 
in the emission (absorption) process of soft Fermi-excitations by a hard particle, this particle 
changes own statistics (unlike the emission or absorption processes of soft Bose-excitations, as 
it is presented in the quantum average (12.141) ). Nevertheless, from the rigorous definition for 
the statistical average (13. 9 p we see that all nonlinear terms in the distribution functions, are 
canceled in exact, and thereby we lead to the simpler expression (13. 3p . It is curious to note 
that quantum average (13. 9p formally coincides with classical one, which follows from (13. 9p in 
the limit 



It is pertinent at this point to note that, generally speaking, analog of quantity is not available in a 
classical plasma, in contrast to the Debye mass mrt- 
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By the next step expression (13. 6p should be compared with the spectral density that follows 
from the fluctuation-dissipation theorem for soft fermionic fluctuations. Unfortunately, as op- 
posed to the bosonic case we could not find somewhere in the literature the FDT relation we 
need. Therefore, two forthcoming sections will be devoted to derivation of the required relation. 



4 Fluctuation-dissipation theorem for source r][ 



a 



In spite of the fact that we are properly interested in the FDT for soft fermionic fluctuations, 
conceptually this theorem can be defined to a limited extent, more easily for fiuctuations of the 
hard Grassmann- valued source r/^. The fiuctuations of sources r/^ and f/^ are thermodynami- 
cally mutually complementary quantities with the fiuctuations of soft fermionic fields ipl^ and 
ipl^. Therefore, knowing the FDT for the first quantities, it is not difficult to restore the FDT 
for the second ones. 

The most rigorous derivation of the fiuctuation-dissipation relations is based on the tech- 
nique of the spectral representation for correlation functions |T^I21]. In the interests of brevity 
we have chosen here a more simple way without using such a representation (see, for exam- 
ple [HE]). First of all we present the total fermion field in the system under study as a sum of 
two parts 

i^a + (4.1) 

where the spinors tpl^ and \E'^ contain the soft and hard momenta, respectively. 

Let us assume that our many-body system is described by the following Hamiltonian H{t) = 
Hq + Hi with the interaction term (in the Schrodinger picture) 

Hi = j rfx[C(x)C,ext(x,t) +<ext(x,t)r)^(x)], (4.2) 

where ipl^ ext(x, t) and ip]^ ext(x, t) are "small" time-dependent Grassmann- valued external fermion 
fields. The color sources rjai'x) and ^^^(x) are expressed in a certain way through the hard part 
(and 'I'Jj) of decomposition f l4.ip . i.e., in other words, through the wave functions of hard 
half-spin particles of the system. The hat above points to the operator nature of the quantity 
under consideration. In what follows, as a simplification in order not to overburden the subse- 
quent formulae, the hat will be systematically omitted. As usual, we suppose that the external 
perturbation is switched off at t = —oo, i.e. 



Hi 



= 0. 

t=—oo 



Furthermore, we introduce into consideration a statistical operator p. This operator obeys 
the quantum Liouville equation 

^h^ = [Ho + Hl,p]. (4.3) 
10 



Here, the bracket symbol signifies the commutator. As the initial condition we make use the 
Gibbs grand canonical ensemble 

= Po = exp{n - P{H^ - ^lN)}, 

where VL is the thermodynamical grand potential, /3 = l/ksT, and /i is a chemical potential 
associated with a baryon number 



N = J ^'(x)7o^Xx) rfx. (4.4) 

Equation (14. 3 p may be solved iteratively, and in a linear approximation in we have 

t 

p{t) ^Po + ^j [HW - t),p,]dt\ (4.5) 

— oo 

where Hl,{t' — t) = Q^Ho{t'-t)jh^i^-iHo{t'-t)/h^ rjj^g expectation value of the operator ^^^(x) is 
defined by the usual formula 

(^;) = MpO- 

Putting the approximation (14.51) into the last expression and making use the cyclic invariance 
of the trace, we have 

{rf^) = «)o + ^y"^(t-t')Sp(po[r/;(x,t),if,!(t')])rft'. 

~oo 

Here, ?7^(x, t) = Q^^ot/h^i^^^-^ ^-iHot/h jg ^j^g operator in the Heisenberg picture with respect 
to the Hamiltonian Hq and (. . .)o = Sp(po . . . ) is an averaging over the Gibbs distribution of a 
global thermal equilibrium configuration. By virtue of local color neutrality and zero fermionic 
number, (?7o)o = 0. 

Further, substituting the external perturbation (14. 2 p into the preceding equation, we find 
iVa) = J Gif\x - x')^ia^') dx' + J 4^a^')F^x' - x) dx\ X = {t, x), 
where we have introduced the double-time Green's functions 

GiTi^ -^') = ^ (^(t - t') Sp(Po{r/L(x, t), r/^(x', t')}) , 

Fl^lW - ^) = ^(^ - t') Sp(Po{r/L(x, t), r^^(x', t')}) . 

Here { , } denotes the anticommutator. The first function represents the usual retarded prop- 
agator for the f] field, the second one is the anomalous propagator. The anomalous Green's 
function in the situation under consideration equals zero. In addition to the retarded Green's 
function we also introduce the double-time advanced Green's function 

GiTix - x') = -^e{t' - t) Sp(po{r]L(x,t),r/^(x',t')}) • 
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To obtain the required fluctuation-dissipation theorem following A.I. Akhiezer and S.V. 
Peletminski [6j, we introduce into consideration a new pair of the Green's functions more 
simple in a structure: 

Jfo.i^-x') = m^',t')vU^,t))o (4.6) 

and 

:7'S(^'-^) = «(x,t)r/^(x',t'))o. (4.7) 
Let us define a relation between these two functions. For (14. 6 p we have 

where in turn 

PoWaix)po = e^(^«-'^^)r7^(x, t) e-^(^«-^^) = e-^'^^r7^(x, t - th/3) e^^^. 

The source is generally an odd function of the operator \E'^ and an even one of the operator 
In our most simple situation, the linearly depends on \I'^, and it does not depend on 
at all. By using the definition of baryon number operator (14.41) and the anticommutation 
relations for the \E'-operators, we have 

e-^'^^r/j,(x, t - ih(5) e^"^ = ?/^(x, t - ihp) e^", 

that immediately results in the following relation between fl4.6p and fl4.7p 

J:^(x' - X, t' - t) = J-j:(x - x', t - t' + zh/3) e-^^, 

or in terms of Fourier components it reads 

j:^(-q,-c.) = J-^:(q,u;)e^('--'^). (4.8) 

The fact of appearance of the chemical potential in an explicit form on the right-hand side of 
(14. 8 p is rather important in subsequent discussioij^. 



^In a number of monographes devoted to quantum statistical mechanics (see, for example, [35]) the gener- 
alized Hamiltonian Hq = Hq — fiN is introduced instead of the Hamiltonian Hq. This leads to the fact that 



on the right-hand side of the quantum Liouville equation (j4.3|) . additional term proportional to /x, appears, i.e. 
instead of (|4.3p now we will have 

^h^ = [Ho+Hlp]-^^[N,p]. 

Repeating the reasoning resulting in (|4.8p . it is not difficult to see that instead of (|4.8p in this case we will 
have ^^^(— q, — w) = J7|^(q, oj) e'^''", i.e., the chemical potential has disappeared. It will be entangled in the 
definitions of the correlators (j4.6l) and (|4.7p in a complicated implicit way. We notice that practically the same 
problem of two possible ways of constructing the fermion propagator at finite temperature and density was 
discussed by A. Niegawa in [23J. 
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The retarded and advanced Green's functions are expressed in terms of the correlation 
functions (14. 6p and (14. 7p as follows: 



gS^^'\x - = ± ^ oi±it - 1')) {J%{x' - ^) + JU^ - x' 

Further, in terms of Fourier components the above expression, in view of Eq. (I4.8p . can be 
written as 



whence, in particular, it follows that 

- G^T^^^^) = i i^'^"^''^ + i)^£(q,^)- (4.9) 

The final step in a derivation of the FDT is introducing into consideration the correlation 
function for fluctuations of the color source r^^. As an initial expression we set the following 
one: 

eU^-x')^^-{K{x),%{x')])o. (4.10) 

In this definition we have considered the fact that (^7^(a;))o = {Va{x))o = 0. The commuta- 
tor under the average sign is reflecting statistics of the ?7^(x) source. We believe that in the 
semiclassical limit the operators r^^ and r^^ turn to classical Grassmann- valued functions. Notice 
that the function (I4.10p possesses the property evident from its definition 

(S^^'(x-x'))^ = 7°S^''(x'-x)7°. 

Furthermore, we present (I4.10p in terms of the correlation functions (14. 6p and ( 14. 7p 

or by using its Fourier components, considering (14. 8 p 

S^,(q,^) = ^(e^^'^-'^^-l)j-£(q,c^). 

After substituting J7^^(q, w) from the above equation into (14. 9p . we arrive at the desired 
fluctuation-dissipation relation 

^%{q,oj) = Uh tanhl/3(;k. - ^^) {g^S'' {c^^u) - G'^^^H^.^)] ■ (4.11) 



We can present the given relation similar to bosonic case ( 12. 9p . Indeed, the partition function 
for bosonic oscillator has the form 



oo 111 

Z, = X:exp{-/3;z..(n + -)} = --^ 



h(/3;za;/2)' 
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and thus its mean energy is 



Qb{oo, T) = — TT-r IhZb = - hio coth j 



dp 2 \2kBT J ' 

as it was defined in f|2.1Up . On the other hand, in the case of fermionic oscillator, at finite 
chemical potential, the partition function has the form 

Zf= exp^-f3{huj - fi)(n - = 2 cosh{f3{nLo - fi)/2) 

n=0,l 

and its mean energy reads 

^, ^, 9,^ 1,, , ^ f hjj — li 

fi, I ) = — T-r In Zp = — [hu — fi) tanh' 



dl3 2' V '^^bT 

Taking into account the preceding expression, the FDT fl4.11l) can be rewritten as follows: 

5 FDT for soft fermionic fluctuations 

After providing a detailed derivation of the FDT for fluctuations of hard source, we move on to 
one of the main goals of this article, namely, to the formulation of the FDT for soft fermionic 
fluctuations. We can obtain the required FDT from the previously derived relation fl4.1ip by 
using the fact of mutuality between the fluctuations of generalized forces (in this case these 
are the sources and fjU and generalized stationary and homogeneous fields in a system (in 
this case these are the soft fermionic fields ipl^ and ipU. This fact was already mentioned at 
the beginning of Section 4. Now we should understand as the perturbation the following 
expression, instead of (14. 2p : 



Hi = j rfx[r/l_,,t(x,t)V^^(x) +V^^(x)r7i^,,t(x,t) 



where f/^ ext(^5 ^) ^"^^ Va cxt(^5 ^) pl^Y ^ ^ole of the generalized Grassmann- valued external forces 
that couple to the operators of soft fermionic fields ipai'x.) and ^/^^(x). For these soft fields we 
introduce the retarded and advanced Green's functions 

Sif'^'^x - x') = ±^0{±{t - t')) Sp(po{C(x,t),^^(x',0}) . (5.1) 

Here, the average is taken over the Gibbs distribution for a system consisting of the medium 
and the radiation of soft fermionic excitations, being at thermal equilibrium with the medium. 
The '^/'-operators are understood as averaged over physically infinitesimal volume with a given 
arrangement of all the hard particles in it. Thereby we restrict ourselves to considering only the 



14 



long- wavelength part of the fermion radiation. This enables us to express the Green's functions 
in terms of macroscopical characteristics of the medium: in this case in terms of the retarded 
and advanced quark propagators in the HTL-approximation [24j. 

The system under consideration can be regarded as involving two subsystems, which by the 
convention can be called the hard and soft ones. As a result of thermal fluctuations there exists 
the fermion number 'pumping over' from one subsystem to another and vice verse. In this 
regard the number of hard particles carrying a half-integer spin (more exactly, the difference of 
the numbers of hard particles and antiparticles) and soft quasiparticles which also obey Fermi 
statistics, is not separately conserved but only its sum is. The chemical potentials associated 
with the hard and soft subsystems subject to thermal equilibrium, are connected by the relation 

/^hard = /^soft = (5.2) 

Let us briefly explain what we understand by mutuality between 77^ and z/'a [13 • Let 
fluctuations in any continuous system be described by some stationary and homogeneous fields 
^^^\x,t), i = 1, . . . , N and let /'^*)(x, t) be a volume density of the generalized force fields, 
conjugate in Lagrange's sense with the generalized coordinates ^'^*^(x, t). The mean power 
dissipated in the volume V under the action of the /^*^(x, t) forces, is written as 

V 

In the paper we have suggested the formula for emitted radiant power of soft fermionic 
excitations. In particular, in the temporal gauge we have 

-r/2 V 

This formula was defined in [25] within semiclassical approximation and therefore the functions 
on the right-hand side are classical Grassmann- valued ones. Comparing (Qf) with (15.31) . we 
can draw the conclusion that if as stationary and homogeneous fields we choose 

e«(x,t) = {^;(x,t),^;(x,t)}, 

then the volume field density of the generalized forces /^*^ here will be 

/«(x,t) = K(x,t),-r/^(x,t)}. 

We believe that such a classification is valid both in (semi) classical and in quantum cases. 
In the previous section for the generalized force, the fluctuation-dissipation relation 

ivW.U = ^ tanh ^-/3{hu - f^){Gif'\ci, uj) - Glf^{c,, u)] (5.4) 
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has been obtained. Here, the (?7a?7^)ajq function represents Fourier component of the correlation 
function for fluctuations of source, Eq. f l4.10p . From the condition of mutuality between the 
generalized forces and coordinates in thermodynamically equilibrium system, a similar relation 
for soft part of the total fermionic field of the system in question, is immediately follows: 

ii^a^U = ^(^^^ tanhl/3(;k^ - f,){sif'\ci,u) - Sif'\ci,co)] , (5.5) 

where now the Green's functions S^^'^^^'' are given by Eq. (15. ip . For global equilibrium system 
the condition of mutuality suggests that these Green's functions are inverse to the functions 
G^'^^^'^ on the right-hand side of ( 15. 4p . 

Further, we will be interested in the semiclassical (long- wavelength) limit of relation (15. 5p . 
In this case the retarded and advanced propagators should be chosen in the HTL-approximation 
according to the rule 

h 

Here, the factor has arisen owing to a choice of our initial definition of the Green's functions 
5'(^,^)_ Within the approximation in question the following useful relation 

*^(^)(g) = /(*5(^)(?))V = -*5(^)(-g) 

is valid. 

Taking into account the expression above in the limit ^ — )■ 0, the FDT (15. 5p results in 

WJiU = (g) + *Si^\-q)}, (5.6) 

where an explicit form of the quark propagator *S^^\q) is given by Eqs. (B.5) - (B.8). Now 
it is necessary to consider the soft spinor fields ipl^ and ■?/'^ are classical Grassmann-valued 
functions. From the relation (15.61) it immediately follows that in order for the spectral density 
for fluctuations of soft fermionic field to be non-zero (in the semiclassical approximation), the 
chemical potential fi of the system under consideration must be non-zerc|f|. 

By using an explicit form of the retarded quark propagator, it is not difficult to obtain 

*S^''\q) + = 22{lm(*A+(g))/i+(q) + Im(*A_(g))/i_(q)}. 

Taking into account the relationship 

Im(*A±(g)) = Im(*Ai^(g))*rA±(g)|2 

and the fact that by virtue of (B.7) and (B.8) one has 

Im(*A^i(g))* = ^(lT^)^^(q^-c.^), 

^ In quantum case (|5.5I) a similar statement takes place for the static limit w — > 0, if the exact propagators 
S" (q, w) have no a singularity on uj in this limit. 
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we obtain the final form of the FDT relation 



ii^>iU = ^ (5.7) 

.2 , ,2\ 



instead of Eq. (15. 6 p . 

Comparing the expression (15. 7p with that for the spectral density (13. 6p obtained within our 
simple model, we see that they differ from each other by the tanh(/3/i/2) factor. Since the 
derivation of expression (15. 7p is more fundamental, it is natural, therefore, to assume that we 
have overlooked something in deriving (13. 6p . Let us return to the procedure of average in (13. 4p . 
Careful analysis of (13. 4 p points out that in this average the following simple fact has not been 
taken into account, namely, the system under consideration is that with a varying number of 
particles and antiparticles. In our case a baryon number N that for massless thermal quarks 
and antiquarks can be presented as 

N = NfNj^^[—, — ^— — ^— 1, (5.8) 

(here, e = |p|), is a measure of change of the particle (and antiparticle) number. The spec- 
tral density (13.61) is generally calculated at some fixed A^. To derive the correlation function 
(V'q(9)'0^(~Q'')) considering a change of A^, we use a trick suggested by R.P. Feynman [26] in 
computing the free energy of the grand canonical ensemble. 

Following pn], we multiply the right-hand side of Eq. (13. 4p by some weighting factor of the 
type: a^Q'^^^ ■ The weighting factor accounts for different probabilities of different numbers 
of particles and antiparticles in the system. However, in our case as distinct from [26], the 
baryon number (15.80 runs not only over all positive values, but also over all negative values by 
virtue of the above-mentioned possibility of exchange of the fermion number between the hard 
and soft fermion subsystems of a quark-gluon plasma. This implies that instead of the factor 
Q-Pi^N need to take e"''^'^'. Further, a careful distinction must be made between even and 
odd number of fermions. As is known (see, for example, [27j) this is achieved by introducing 
the fermion number operator (—1)^ (the factor (—1)^ in our case), which is identified with the 
qn coefficient. Taking into account all the above-mentioned, instead of (13. 4p we have to write 
now 



27r2 



X / ^ CSiq)x)^ ix *Si~q')), 5(a; - V q) 5(g - q'). 



Considering 



+ oo + oo ^ 

J2 (_l)^e-^^l^l = 1 + 2 ^(-l)^e-^^^ = tanh(^/3/i), (5.9) 



iV=-oo N=l 
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we conclude that to account for the fact of a change of the particle number in spectral density 
( 13. 6p . it is sufficient to multiply it by the factor tanh(/3/i/2). This results in an identical 
coincidence between the spectral densities fl3.6l) and (15.71) . The only difference is that for the 
first spectral density, the sum (15. 9p converges only for /i > 0, while (15. 7p is true for any values 
of the chemical potential. 

In the case of < as the weighting factor it is necessary to take e"'"^'^'^', where now we 
have to set = —(—1)^. Then instead of the sum (15. 9p we will have 

+ oo + oo ^ 

_ ^ (_i)^e^^l^l = -1 - 2 ^(-l)^e^^^ = tanh(- 

N=-oo N=l 

and in doing so we return again to the same factor, but with the negative chemical potential. 

6 Conclusion 

In this paper we have proved two independent derivations of the spectral density for soft 
fermionic excitations in an equilibrium quark-gluon plasma. Careful consideration of all features 
that are common to excitations (both soft and hard) obeying Fermi-Dirac statistics, has allowed 
us to obtain total coincidence of the final expressions. This in itself is a remarkable fact, since the 
initial premises of these two approaches, are supported on different physical grounds. Another 
very intriguing feature is rather strong formal similarity of our first derivation with the bosonic 
case. The latter was considered in Section 2 for comparison. 

The result obtained in this work appears to be sufficient justification for our pseudoclassical 
approach jUIS] in a description of the dynamics of hard and soft excitations of a hot QCD 
plasma obeying various statistics, at least on the level of thermal fiuctuations (at soft scale 
gT) , within perturbation theory. The next step should be an extension of the model suggested 
in [1112] to a description of highly excited states, i.e. QGP far from thermal equilibrium and 
beyond the limits of perturbation theory (strongly coupled plasma). 
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Appendix A 



In Section 2 we make use the following formula of the integration with respect to commuting 
color charge Q": 

jdQQ^Q'= C = Q,Q,G, (A.l) 

where as the integration measure dQ the following expression is meant [T ^ [T ^ [25]: 

dQ = d^'^Q cr 6{Q''Q'' - ga) • • • • (A.2) 

The constant cr is fixed by the normalization condition J dQ = 1. With the constant q2 one 
relates the so-called quadratic Casimirs. The dots on the right-hand side of (A.2) denotes delta- 
functions ensuring the conservation of the higher Casimirs for the color groups SU (Nc), Nc > 3. 

By the quadratic Casimirs in quantum field theory, the numerical constants of the unit 
operator / of proper dimension on the right-hand side of the expressions 



= CfI, 
T^T" = CaI, 



are commonly taken, where (T"')^^ = —if"''"^, Cp = {N^ — l)/2A^c; and Ca = N^. To provide 
the correct expressions in integrating with respect to the color charges as the quark and gluon 
quadratic Casimirs in (A.l), we will mean 

= tr(t"t") = TpdA, (A.3) 

Cf ^ = tr(T'^T'^) = CacIa. (A.4) 

Here, dA = — 1 is dimension of the gauge group and Tp = ^ is the index of the fundamental 
representation. For such a choice of the constants Cg^'', in particular, we derive the correct 
expression for the Debye screening mass (I2.15p . The constant q2 in measure (A.2) is simply 
C2^K This is exactly what is given in the paper |28] . 



Appendix B 

In this Appendix all necessary formulae for the gluon and quark propagators in the hard 
thermal loop approximation and also some formulae for the integration over solid angle we use 
throughout this work, are given. 

The HTL-resummed gluon propagator in the covariant gauge is defined by the following 
expression: 

*V^,ik) = -P^M *^\k) - Q^,{k) *A\k) + ^D^,{k)A%k), (B.l) 
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where A^{k) = l/k^; *A*''(A;) = 1/(A;2 - n*''(A;)), U^{k) = lU^"'{k)P^^{k), and n'(A;) 
U.'^'^(k)Qfj,,y{k) with the gluon self-energy n^j,(/c) 



(B.2) 



^ is a gauge- fixing parameter. An exphcit form of the Lorentz matrices P^y{k), Qixu{k), and 
D^y^k) in the chosen gauge, is 



Pfj-fik) — g/xf D^vik) QiJ,v{k\ Qixvik) — -211 \ ' ■^iJ.i^i^) — 



(B.3) 



Uf,{k) = k'^u^ - k^{k ■ u), 

where Ufj, is the global four-velocity of a non-Abelian plasma. The Lorentz matrices possess 
evident properties of ordinary projectors 

P2 q2 ^ g J^2 ^ pQ^ pjJ^QJJ^ 0. (B.4) 

Further, the medium modified quark propagator *S{q) has the following form: 

*S{q) = /i+(q) *AM + /i_(q) *A_(g), (B.5) 

where the matrix functions 

/i±(q) = (7°Tq-7)/2 (B.6) 
with q = q/|q|, are the spinor projectors onto eigenstates of helicity, and 



1 



uj^[\q\ + SE±{q)] 
are the 'scalar' quark propagators, where in turn 



Vlq| 



(B.7) 



(B.8) 



with 



In 



1 + z 



1 - z 



in9{l - \z\) 



are the scalar quark self-energies for normal (+) and plasmino (— ) modes. For the h±{q) 
matrices the following identities 



^±(q)^±(q) = 0, /i±(q)/iT(q)/i±(q) = /i±(q) 

are valid. 

Finally, we give an explicit form of some simple integrals over solid angle dflv'- 



(B.9) 



r 

J 1^ 



(B.IO) 



An 

where /c* = /cVlkl. 



v'^v^5{uj — V • k) 
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Appendix C 



In this Appendix we discuss a problem of construction of the integration measure with respect to 
anticommuting color charges ^* and 6'^*. These charges obey the following dynamical equations 

^ + ^gv^A;{t,^rt){tTO^{t) = 0, 9}, = 9\t)\ 

de%t) ^^-'^ 

—±l-^gv^A^^{t,vt)e^\t){ty' = 0, 0l' = 0^\t) , 
at ^ t=o 

where v'^ = (1, v). Here, we have neglected the interaction of the Grassmann charges with soft 

ip-fields induced by thermal fluctuations in a hot non-Abelian plasma. This system of equations 

permits an quadratic integral of motion 

e^'it)9\t) = 9l'ei,^cP, c = Q,Q,G. 



The integral of motion is analog of the invariant Q°'(t)Q°-(t) = QqQq ( = C'g'''') given in Appendix 
A and it as well as the invariant Cg''^ depends on the type of hard parton to which the Grassmann 
color charges concern. 

As a definition of the Grassmann color charge integration measure we put the following 
expression: 

dOde^ = (j\d9'd9^'^f{9^9), (C.2) 

where f{9^9) is unknown for the time being function of the invariant 9'^ 9 = 9'^^9^ . This measure 
is invariant with respect to the involution f. The function / must formally play a role of a 
Dirac (5-function as it takes place, for example, in the definition of measure for the usual color 
charges, Eq.(A.2), i.e., 

f{9^9)^5{9^9-Cf^). 

Let us look for an explicit form of this function in a way of a (finite) expansion in powers of 9'^ 9 

fio^o) = h + h e^e + h{9^9f + ... + h^.,i9^9f^-' + fNA0^9f^. (c.3) 

Here, /„, n = 0, 1, . . . , A'^c are some constants depending on the group invariants. Below we use 
the standard rules for the integration over the Grassmann variables 



d9' = j d9^' = 0, 

Let us require that for each function of the form g = g{9^9) the relation 



jd9d9^g{9^9) = g{&P). (C.4) 
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was fulfilled. The function g{d'^d) represents a finite polynomial in O'^O similar to (C.3). Making 
use the integration rules over the Grassmann variables, it is not difficult to obtain that the 
relation (C.4) will be fulfilled if the function f{d^d) has the form 

f[e^e) = ^l{cfY^+{cfY^-'e^e+ ... +cf\9^9)''^-^ + {e^ef^. (c.5) 

As the simplest consequences of this fact we have the normalization condition 

j dOdO^ = 1. 

and also the integral 

j d9d9^ 9^'9^ = j S''- (C.6) 

It is worth noting here that the integral (C.6) enables us to calculate the other more complicated 
expressions involving the anticommuting color charges of different particles. So when doing 
the calculation of the probability of soft quark bremsstrahlung j2], we are faced with a color 
structure of the following type: 

Cif^ ^ i9\t'^92)i9lt''0i)- 

Here, the labels 1 and 2 concern to the first and second hard particles involved in the scattering 
process in question. The integration over the Grassmann color charges of these two hard partons 
of this color factor is easily carried out by the general formula (C.6) 



d9id9l d92d9lclf^ = -tr(n") I ^ 1 I ) = -C. 



/r^(2)\ /^(i)^(2)- 



(1-2) 

The explicit value obtained for the color factor C^g , in exact coincides with a similar one in 
the paper [2j derived from a basically different reasoning by means of much more cumbersome 
calculations. 

It is also easy to calculate more nontrivial expression of the type jd9d9'^ 9'^'^9^9'^^9K Obvi- 
ously, the given expression has to be antisymmetric with respect to a permutation of indices 
i ^ k and j ^ L Therefore, we have 

Jd9d9^ 9^'9^9^^9^ = a{6'^^^ - 5"^ 5^^), (C.7) 

where a is some constant. A rather easy calculation leads to 

" N.xN,-iy 

The integral 

jd9d9\9h'^9){9H'9) = [tr^r - tr ^t^] = -^^^^^ 
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is a particular consequence of relation (C.7). In a similar way we can calculate the other more 
complicated integrals. 

Thus, setting into consideration the integration over the Grassmann color charges enables us 
to automatize fully the procedure of calculation of various color factors which appear when doing 
the determination of the scattering probabilities involving hard and soft fermionic excitations. 
In doing so, we practically achieve here, perfect analogy with the purely bosonic case [30j. 
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